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ABSTRACT 

We use high-resolution UV spectra of the radio-quiet QSO H1821+643 (zcm = 0.297), 
obtained with the Space Telescope Imaging Spectrograph (STIS) and the Far Ultraviolet 
Spectroscopic Explorer (FUSE), to study the ionization and metallicity of an intervening O VI 
absorption line system at Zabs — 0.1212. This absorber has the following notable properties: 
(1) Several galaxies are close to the sight line at the absorber redshift, including an actively 
star-forming galaxy at a projected distance of 144 h^^ kpc. (2) There is a complex cluster 
of H I Lya absorption lines near the O VI redshift, including at least five components spread 
over a velocity range of ~700 km s^^. (3) The strongest Lya line in the cluster appears to be 
composed of a mildly saturated component with a typical fo— value blended with a remarkably 
broad component with & « 85 km s~^. (4) The O VI absorption is not aligned with the strongest 
(saturated) H I absorption, but instead is well- aligned with the very broad component. (5) The 
only detected species (at the 4a level) are O VI and H I despite coverage of strong transitions 
of abundant elements (e.g., C II, C III, and C IV). Based on these constraints, we find that the 
absorption line properties can be produced in coUisionally ionized gas with lO^ '^ < T < 10^-^ K 
and —1.8 < [0/H] < —0.6. However, we find that photoionization is also viable if the pathlength 
I through the absorbing gas is long enough; simple photoionization models require 85 < / < 
1900 kpc and —1.1 < [0/H] < —0.3. We briefly discuss how observations of X-ray absorption 
lines due to O VII and O VIII could be used, in principle, to break the ionization mechanism 
degeneracy, and we conclude with some comments regarding the nature of O VI absorbers. 
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1. Introduction 

Any successful cosmological model must account for the quantity, distribution, and physical state of 
the baryons in the Universe. While this is likely a more tractable challenge than understanding the more 
exotic non-baryonic dark matter and dark energy (i.e., the cosmological constant), there is a "missing 
baryon problem" which has been known for some time. In brief, inventories of the well-observed baryons 
in the nearby Universe (including stars, cool neutral gas, and X-ray emitting galaxy cluster gas) find 
total baryon densities, in units of the closure density, of fib = p/pc ~ 0.006 (assuming Hq = 75/175 km 
s~^ Mpc~^, see, e.g., Persic & Salucci 1992; Fukugita et al. 1998). This is far less than the expected 
value based on Big Bang nucleosynthesis and deuterium measurements (Buries & Tytler 1998): f^b/if 00 — 
0.019±0.001 ( = 0.034±0.002 for hioo = 0.75). Ironically, the baryons have been more readily accounted 
for in the distant Universe than in the low— z Universe: at high rcdshifts, the "Lya forest" - photoionizcd, 
diffuse intergalactic gas at T ~ 10'' K - appears to contain the vast majority of the baryons predicted by 
the deuterium observations (Ranch et al. 1997; Weinberg et al. 1997). At low redshifts, the cool Lya 
forest is less dominant but still makes a substantial contribution to the baryon inventory. However, even 
including recent estimates of the baryonic content of the low— z Lya forest (ShuU, Stockc, & Pcnton 1996; 
Penton, ShuU, & Stocke 2000), the census of baryons at the present epoch falls well short of the expected 
fib- Consequently, cosmologists are confronted with an important question: where are the missing baryons 
at the present epoch? 

A promising possible answer to this question has been provided by cosmological simulations of the 

growth of large scale structure. According to hydrodynamic simulations of structure formation, at the 
present epoch 30—50% of the baryons (by mass) are located in low-density, intergalactic gas which has been 
shock-heated to 10^ - lO'^ K (Cen & Ostriker 1999a; Dave et al. 1999, 2001). Following Cen & Ostriker, 
we refer to this 10^ — 10'' K gas as the warm/hot intergalactic medium to distinguish it from the hotter 
gas observed in rich galaxy clusters. Such gas is challenging to observe. Its X-ray emission is extremely 
difficult to detect, especially at lower energies where confusion due to foreground emission and absorption is 
substantial. Current X-ray observations do not exclude this hot IGM baryon reservoir (Kuntz, Snowden, & 
Mushotzky 2001; Phillips, Ostriker, & Cen 2001), and there are some indications of large-scale filaments of 
diffuse gas delineated by X-ray emission (e.g., Scharf et al. 2000; Rines et al. 2001). However, Voit, Evrard, 
& Bryan (2001) and Croft et al. (2001) have demonstrated that it is difficult to distinguish between X-ray 
emission from diffuse intergalactic gas outside of virialized structures (as predicted by the cosmological 
simulations) and X-ray emission from virialized groups and clusters. The gas may be a significant baryon 
repository in either case, but the distinction has important physical implications, e.g., regarding the relative 
importance of gravitational vs. nongravitational heating processes. 

QSO absorption lines provide an alternative means to test the prediction that a significant quantity 
of baryons are in the warm/hot IGM. Absorption lines such as the O VI AA1031.9,1037.6 and Ne VIII AA 
770.4,780.3 doublets provide sensitive probes of 10^ — 10^ K gas in coUisional ionization equilibrium (Verner, 
Tytler, & Barthcl 1994). Searches for the O VI absorption lines at Zabs^ 0.3 with the Space Telescope 
Imaging Spectrograph (STIS) on board the Hubble Space Telescope (HST) have revealed that the number 
of intervening O VI absorbers per unit redshift (dN/dz) is remarkably high (Tripp, Savage, & Jenkins 2000; 
Tripp & Savage 2000). Furthermore, with reasonable assumptions about the metallicity and ionization of 
the gas, these studies have shown that O VI absorption systems are probably an important baryon reservoir. 

Since the O VI ion fraction peaks at T ~ 3 x 10^ K in collisional ionization equilibrium, these O VI 
absorption studies provide tantalizing evidence that the warm/hot IGM is an important baryon repository 
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at 2; ~ 0. However, there are several issues which must be addressed in order to build a compelling case. 
First, the number of O VI systems (and sight lines) in the published STIS observations is rather small, 
and consequently quantities such as dN/dz and the cosmological mass density have large uncertainties. 
More observations are needed, and new programs to build the sample are underway. Second, the O VI 
lines could arise in gas which is not coUisionally ionized, but rather is photoionizcd by the UV background 
from QSOs and active galaxies. Or, the gas could be coUisionally ionized but out of equilibrium, e.g., 
because it is able to cool more rapidly than it can recombine. Photo- or non-equilibrium ionization is 
favored in some of the O VI systems because they have narrow associated H I lines (e.g., Tripp & Savage 
2000; Savage et al. 2001) which would not arise in equilibrium coUisional ionization. Alternatively, these 
could be multiphase absorbers. Of course, in some cases photoionization and coUisional ionization could 
be comparably important. If the gas is photoionized, then there are concerns about double-counting in 
the baryon inventory (see below). Third, the metallicity of the absorbers is poorly constrained. If these 
O VI systems originate in pockets of high metallicity gas, then there is less hydrogen associated with the 
absorbers and their baryonic content is lower [f2(,(0 VI) oc (0/H)~^]. Finally, the environment in which 
these absorbers are located must be scrutinized. Are these absorbers found in unvirialized galaxy filaments 
or even galaxy voids, or do they arise in virialized structures? Are these systems due to gas which is 
gravitationally shock-heated when gas accretes onto large-scale structures, or could they be heated by 
nongravitational processes such as supernova-driven winds? 

In this paper we are primarily interested in the ionization and metallicity of the gas. Specifically, we 
use high-resolution UV spectra obtained with STIS and the Far Ultraviolet Spectroscopic Explorer (FUSE) 
to study the ionization and nature of a particular absorption system at Zabs = 0.1212 in the spectrum of 
the radio-quiet QSO H1821-I-643 (2;em = 0.297). While some O VI absorbers are apparently "intrinsic" (i.e., 
close to the QSO itself, see Hamann & Ferland 1999), the 111821+643 absorber at 2;abs = 0.1212 is certainly 
an intervening system: it is highly displaced from the QSO redshift (At; w 43,400 km s~^) and is located 
at a redshift where there are several galaxies close to the line of sight, including a luminous galaxy at a 
projected distance of 144 h^^ kpc which shows [O II], [O III], and H/? emission lines indicative of active 
star formation (Tripp, Lu, & Savage 1998; Bowen, Pettini, & Boyle 1998). The ionization mechanism is a 
crucial issue because if the O VI absorbers are photoionized and cool rather than coUisionally ionized and 
hot, then they still represent a substantial baryon reservoir (see Tripp & Savage 2000), but in this case they 
may reveal the same baryons counted in Lya studies such as Penton et al. (2000) [note that all intervening 
O VI systems that we have detected so far are also detected in the H I Lya transition]. Apart from the 
baryon census, the metallicity and physical conditions of the gas are interesting in their own right. 

The paper is organized as follows. The STIS and FUSE observations and data reduction are described 

in §2. Section 3 summarizes the absorption line measurement techniques. Constraints are placed on the 
physical conditions and metallicity of the gas in §4 including a comparison of the absorber properties 
required in coUisionally ionized and photoionized scenarios and some comments on X-ray absorption lines 
of O VII and O VIII. We discuss the results and summarize our conclusions in §5. Throughout this paper 
we assume i/o = 75 km s^^ Mpc^^ and = 0.0 imless otherwise indicated. We also report heliocentric 
wavelengths and redshifts; in this direction, wlsr = Vh.e\io + 16 km s~^ assuming the standard definition of 
the local standard of rest (Kerr & Lynden-Bell 1986). 
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2. Observations and Data Reduction 

This paper makes use of observations of H1821+643 taken with STIS in both the echelle and first-order 
grating modes, as well as shorter-wavelength observations made with the FUSE satellite. In this section we 
briefly describe the observations and data reduction; see Woodgate et al. (1998) and Kimble et al. (1998) 
for information on the design and performance of STIS, and Moos et al. (2000) and Sahnow et al. (2000) 
regarding the FUSE design and performance. 

2.1. STIS Spectroscopy 

H1821-h643 was observed with the medium resolution FUV echelle mode (E140M) of STIS on 1999 
June 25 and 2000 March 31 for a total integration time of 50.93 ksec. This STIS mode uses the FUV 
MAMA detector, which generally has very low dark background counts. While this detector also has an 
amoebic region that shows elevated dark counts when the detector temperature increases (see Figure 1 
in Brown et al. 2000), the dark count rate was small compared to the source count rate for all of the 
observations employed here. The observations made use of the 0'.'2 x 0'.'06 slit resulting in a line spread 
function with minimal wings (see Figure 13.87 in the STIS Instrument Handbook, v4.1). This STIS echelle 
mode/slit combination provides a resolution oi R = A/AA ~ 46,000 (FWHM sa 7 km s^^, see Kimble et 
al. 1998) and wavelength coverage from ^1150 to 1710 A with four small gaps between orders at A > 1630 
A. Several STIS first-order grating observations of 111821-1-643 were also obtained on 1999 June 24. Most 
relevant to this paper is the STIS G230M observation extending from 1724 to 1814 A, which covers the 
C IV doublet at Zabs = 0.1212. This first-order grating spectrum has a resolution of ^30 km s^^ (FWHM). 

The data were reduced with the STIS Team version of CALSTIS at the Goddard Space Flight Center. 
The individual echelle spectra were flatfielded, extracted, and wavelength and flux calibrated with the 
standard techniques. Then a correction for scattered light[| was applied using the method developed by 
Bowers et al. (2001), and the individual spectra were combined weighted by their inverse variances averaged 
over a large, high S/N region. Finally, overlapping regions of adjacent orders were also coadded weighted 
inversely by their variances.^ The first and last ten pixels in each order were not coadded to avoid spurious 
pixel values often seen near order edges. 

The STIS G230M first-order grating observation covering the C IV AA1548. 2, 1550.8 doublet at Zabs 
= 0.1212 was flatfielded and wavelength calibrated in the normal way, but the spectrum was extracted 
with the optimal method of Robertson (1986). Two 50-pixel wide regions centered 52 pixels away from 
the spectrum (one on each side) were used to determine the background. The two G230M spectra were 
also combined with weighting based on signal- to-noise. To ensure that the wavelength scale of the G230M 
spectrum is aligned with that of the E140M echelle spectrum, the Milky Way Ni II A1741.5 and A1751.9 
lines covered in the G230M spectrum were compared to the Milky Way Ni II A1317.2 and A1370.1 lines 



^Inspection of the cores of strongly saturated lines indicates that while the scattered light correction is usually quite effective, 
occasionally errors in the flux zero point at the level of a few percent of the continuum are evident. Consequently, we have 
included a term due to the flux zero point uncertainty in the overall uncertainties in the STIS measurements (see §3), and we 
have allowed the profile-fitting code to adjust the zero point by a few percent (for the STIS data only) as a free parameter. 

*For the overlapping regions of adjacent orders, the weighting was determined on a pixel-by-pixel basis. However, a five- 
pixel boxcar smoothing was applied to the error vectors for the determination of the weights so that pixels with larger noise 
fluctuations are not inappropriately over- or under-weighted. 
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recorded in the E140M spectrum. 

2.2. FUSE Spectroscopy 

H1821+643 was initially observed with FUSE on 1999 October 10 and 13; the total integration time 
was 48.8 ksec. Results from this observation have been presented by Oegerle et al. (2000). Subsequently, 
the QSO was re-observed with FUSE on 2000 July 24 with a total integration time of 62.7 ksec. Both 
observations used the large (30" x 30") LWRS aperture. FUSE has four co-aligned telescopes and Rowland 
spectrographs which record spectra on two microchannel plate detectors. Two of the optical channels have 
LiF coatings to cover the 1000-1187 A range, and the other two channels have SiC coatings to cover 
905—1105 A. During the 1999 October observation, only the two LiF channels were successfully aligned, 
but high-quality spectra were obtained with all four channels in 2000 July. By design, most wavelengths in 
the FUSE bandpass are recorded by at least two channels (see Figure 3 in Sahnow et al. 2000). 

The initial (1999) observations were reduced as described in Oegerle et al. (2000). The later 
observations were reduced in a similar fashion with CALFUSE version 1.7.6. The data were recorded in 
time-tag mode and were processed to correct for Doppler shifts due to orbital motion, subtract backgrounds, 
and wavelength and flux calibrate the spectra. Flatfield and astigmatism corrections were not sufficiently 
tested at the time of processing and were not applied. CALFUSE assumes that the background is uniform 
across a given detector, which is an adequate approximation for an object as bright as H1821+643. Oegerle 
et al. (2000) report that the resolution of the 1999 October spectrum is 20—25 km s~^. The resolution of 
later observations are somewhat better; the FUSE PI Team finds that the current resolution ranges from 
~17 to 25 km s~^ (e.g.. Savage et al. 2001). In both cases the resolution depends on the wavelength and 
channel used to record the data. Because of the varying resolution of data from different channels and 
observation dates, there are differing views on how to best make use of multiple detections of a particular 
line. We discuss our strategy regarding this issue in §3. 

While the wavelength scale of the 1999 October data can show large errors (10-30 km s~^) over small 
intervals, this problem was subsequently corrected and is not evident in the 2000 July data. The relative 
dispersion solution of the 2000 July data is accurate to ~6 km s~^. However, there can still be an error in 
the zero point of the wavelength scale, e. g., due to imperfect centering of the target in the aperture. To set 
the zero point of the FUSE data, the multicomponent Milky Way Fe II A1144.9 line was compared to the 
Galactic Fe II 1608.5 line recorded in the STIS E140M spectrum. Similarly, Lyman series lines in the FUSE 
spectrum at Zabs = 0.225, which show distinctive component structure, were compared to the Ly7 line at 
this redshift in the STIS spectrum. As in Oegerle et al. (2000), the wavelength scale of the 1999 October 
data was corrected on a line-by-line basis by comparison to appropriate analogous lines in the STIS E140M 
observation. 



3. Absorption Line Measurements 

The absorption profiles of detected lines at ^abs = 0.1212 as well as regions of undetected species of 
interest are shown in Figure 0. Standard processing of FUSE data provides highly oversampled spectra, 
so we have rebinned the FUSE data to « 7 km s~^ pixels. The STIS data are optimally sampled, so 
no rebinning has been applied. Remarkably, despite the close proximity of several galaxies (see §1) and 
a relatively strong Lya line, only Lya, Ly/3, and the O VI A1032 transitions are clearly detected at the 
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Fig. 1. — Continuum-normalized absorption profiles of the H I Lya, Ly/3, and O VI A 1032 lines detected at z^bs = 0.1212 
in the spectrum of H1821-I-643 (lower three profiles); the profiles are plotted vs. rest-frame velocity where ti = at z^bs = 
0.1212. The C II A1334.5, C III A977.0, C IV A1548.2, and N V A1238.8 lines are not significantly detected at this redshift, 
and the spectral regions of these undetected lines are shown in the upper portion of the figure. The profile fit which yields the 
component parameters in Table ^ is overplotted on the Lya and Ly/3 profiles with a solid line, and the broad H I component 
discussed in §4 is shown with a dashed line. For reference, the velocity centroid of this broad component is indicated with a 
tick mark above each profile. Note that there are several unrelated lines near the C III and Ly/3 profiles. 
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4ct level or better at this redshift (we argue below that there is evidence of O VI A1038 absorption as 
well, although the transition is blended with an unrelated line). Table |l| lists rest-frame equivalent widths 
(Wr) of the detected absorption lines at ^abs = 0.1212, measured using the methods of Sembach & Savage 
(1992), as well as 4(t upper limits on undetected lines. Note that contributions due to the uncertainty in 
the curvature and height of the continuum as well as a 2% uncertainty in the flux zero point are included 
with the statistical uncertainties in the overall errors in Table |l|. Integrated apparent column densities 
(Savage & Sembach 1991), which similarly include continuum placement and zero point uncertainties, are 
also provided in Table |l|. 

The quantities in Table |l] are integrated over the full velocity range of the main Lya feature but do 
not include the three high- velocity weak Lya lines evident in Figure |l| at w = —425, —179, and 4-270 km 
s^ . However, it is obvious from Figure |l| that even the main Lya feature is composed of at least two 
substantially blended components. To deblend these components and measure line widths, we have used the 
Voigt profile-fitting software of Fitzpatrick & Spitzer (1997) to simultaneously fit the Lya and Ly/3 profiles. 
For the STIS echelle data (Lya), we have accounted for instrumental broadening using the line spread 
functions from the STIS handbook; for FUSE data (Ly/3), we have assumed that the LSF is adequately 
described by a Gaussian with FWHM w 20 km s~^. The profile- fitting results are summarized in Table 
The final model profile is overplotted on the Lya and Ly/3 profiles in Figure |l| with a solid line. 

Two comments on the profile-fitting results are in order. First, the Ly/3 profile shows excess absorption 
compared to the model profile over several pixels at w « 90 km s^^. This is due to noise. Since the Lya 
/—value is a factor of 5 greater than the Ly/3 /—value, and since the equivalent width goes as /A^, this 
"component" would be very prominent in the Lya profile if it were real. We have attempted to find a set 
of components which produce this extra Ly/3 absorption without violating the constraints set by the Lya 
profile, and we have been unsuccessful. Since FUSE is known to have complex and substantial fixed-pattern 
noise (Sahnow et al. 2000), it is not too surprising to occasionally encounter a noise feature of this sort. 
Second, the model profile fitted to the H I lines is not unique. The smoothness of the Lya profile at 
100 <, V <^ 200 km s^^ provides some support for fitting this wing with a single broad component, but a 
comparably good fit can be obtained with several narrower components spread over this velocity range. 

In the case of the FUSE data, we are primarily interested in total equivalent widths and column 
densities integrated across the full velocity range of the main absorption feature (i.e., —110 < v < +210 km 
s~^ in Figure Given the breadth of the absorption features of interest and the moderate differences in 
resolution of the various data, it is reasonable to coadd the various spectra to increase the S/N. Nevertheless, 
to be conservative we have measured integrated equivalent widths and apparent column densities from 
the individual spectra separately, then we have determined the mean of the individual measurements, 
each weighted inversely by its variance. These weighted means are the final integrated quantities reported 
in Table |l|. For profile fitting, on the other hand, we fitted the STIS Lya and FUSE Ly/3 profiles 
simultaneously. Consequently, to make full use of the FUSE observations, we first coadded Ly/3 profiles 
with the usual weighting based on the inverse variances of the individual spectra. We similarly coadded the 
FUSE O VI data for the purposes of the next paragraph. Note that Figure |l] shows these coadded Ly/3 and 
O VI profiles. Since the Ly/3 lines are detected at high significance in the individual spectra, we elected 
to coadd only the Ly/3 profiles from the two LiF channels obtained in 2000 July and thereby preserve the 
resolution as much as possible. However, the O VI lines are detected at lower significance in the individual 
channels, and one of the LiF channels from 2000 July is not useful in the O VI region due to the vertical 
stray light stripe caused by scattered terrestrial Lya emission during spacecraft orbital daytime (see §2.3 
in Sahnow et al. 2000). Consequently, for the O VI lines we have coadded both of the LiF channels from 
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Table 1. Equivalent Widths and Integrated Column Densities of the O VI Absorber at z = 0.1212 



Species 


Ao" 


r 




log iVa^ 


Spectrograph 




(A) 




(ml) 






HI 


1215.67 


0.416 


595±12 


>14.38'^ 


STIS E140M 


HI 


1025.72 


0.0791 


158±14 


14.43±0.04 


FUSE<^ 


VI 


1031.93 


0.133 


90±17 


14.02±0.07 


FUSE* 


N V 


1238.82 


0.156 


<56s 


<13.42s 


STIS E140M 


CIV 


1548.20 


0.191 


<144S 


<13.55s 


STIS G230M 


cm 


977.02 


0.759 


<104s 


<13.2is 


FUSE'' 


CII 


1334.53 


0.127 


<62g 


<13.49s 


STIS E140M 


Si IV 


1393.76 


0.514 


<55s 


<12.79s 


STIS E140M 


Si III 


1206.50 


1.669 


<70s 


<12.5is 


STIS E140M 



''Rest frame vacuum wavelength and oscillator strength from Morton 
(2001) or Morton (1991). 

''Rest frame equivalent width integrated from —110 to 210 km s~^ 

where v = km at Zabs = 0.1212. 

■^Apparent column density integrated from —110 to 210 km s~^. For 
undetected lines, Aa upper limits are derived from the upper limits on 
Wr using the linear portion of the curve of growth. 

''Saturated absorption line. 

''Weighted mean of individual measurements from the observations 
with the LiFlb and LiF2a detectors on 1999 October 17 and 2000 July 
24. 

^ Weighted mean of individual measurements from the observations 
with the LiFlb and LiF2a detectors on 1999 October 17 and the LiF2a 
detector on 2000 July 24. The LiFlb spectrum from 2000 July 24 was 
adversely affected by the vertical stray light "stripe" (Sahnow et al. 
2000) and is too noisy to provide a useful measurement. 

^4(7 upper limit. 

''Based on the LiF2a spectrum from 2000 July 24 only. 
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Table 2. Component Parameters of the Lya and Ly/3 Absorption Lines at z = 0.1212'' 



Redshift*^ v'= b (km s"^) log N 



0.11961 


-425 


35±9 


13.14±0.04 


0.12053 


-179 


23+1' 


12.71±0.08 


0.12112 


-21 


26+1' 


13.93±0.37 


0.12125 


13 




14.04±0.36 


0.12147 


72 


&5+_ll 


13.78±0.17 


0.12221 


270 


38±9 


13.22±0.04 



''Based on simultaneous profile fitting of 
the Lya and Ly/3 lines with the software of 
Fitzpatrick & Spitzer (1997) and the STIS 
E140M Une spread functions from the STIS 
Instrument Handbook. The FUSE LSF was 
taken to be a Gaussian with FWHM w 20 km 
s-i. 

''Heliocentric redshift; the conversion from 
a heliocentric velocity scale to the Local 
Standard of Rest velocity scale is given by 
v^LSR = Vh.e\io + 16 km s~^ assuming the Sun 
is moving in the direction I = 56°, 6 = 23° at 
19.5 km s-i (Kerr & Lynden-Bell 1986). 

■^Component velocity centroid in the rest 
frame of the absorption system where v = 
km s~^ at Zabs = 0.1212. 
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the 1999 observations with the one useful LiF channel (LiF2) from the 2000 observations. We note hat the 
measurements from the individual channels are consistent within their la uncertainties. 

Since O VI has a resonance-line doublet, its absorption lines can usually be unambiguously identified. 
Unfortunately, the weaker 1038 A line of the O VI doublet at z^hs = 0.1212 in the spectrum of H182 1+643 
is strongly blended with the Lyd line from the absorption system at Zabs = 0.225, which may lead some 
readers to question the reliability of the identification of the stronger O VI A 1032 line shown in Figure |^. 
However, since the hy6 line is not strongly saturated, we should be able to demonstrate that the 1038 A 
line is present if the S/N is sufficient, and indeed we do find evidence of the 1038 A line at the expected 
wavelength and strength. Figure ^ shows this evidence. This figure shows the continuum-normalized, 
coadded O VI profiles. A single-component fit to the O VI A1032 line only is shown with a dashed line 
overplotted on the 1032 A profile (upper histogram). The velocity, column density, and Doppler parameter 
from this fit predict a 1038 A absorption line shown with a dashed line shown overplotted on the lower 
histogram in Figure |^. Absorption is clearly present in the lower profile at the expected velocity with the 
expected strength. This gives us confidence that the O VI identification is correct. 



4. Physical Conditions and Metallicity 

In this section we explore the implications of the absorption-line measurements regarding the ionization 
and metallicity of the gas. We first consider coUisional ionization, and then we examine photoionization. 
Before assessing the ionization mechanism, we comment on several interesting features of this absorber: 

1. The O VI absorption occurs at a redshift where there is a cluster of Lya absorption lines; at least five 
Lya lines are present within a velocity interval of ~700 km s^^ (see Figures |l| and|3|). 

2. Furthermore, there are several galaxies within ±500 km s~^ of this redshift with projected distances 
from the sight line ranging from 144 hj^ kpc to 3.1 Mpc (see Tripp et al. 1998). Figure || shows the 
velocities and projected distances of these galaxies with respect to the Lya absorption lines [Tripp et 
al. (1998) estimate that the galaxy redshift uncertainties are ~ 50 — 100 km s~^]. As noted in §1, the 
closest galaxy has an emission line spectrum which suggests that it is actively forming stars. 

3. In addition to an ordinary, mildly saturated narrow line, the main Lya profile appears to contain a 
relatively broad component at v — 72 km s^^ with 6 w 85 km s~^ (see Table 2; this broad component 
is shown with a heavy dashed line in Figure This implies an upper limit on the temperature of 
the gas at this v: T < mb^/2k = 4.3 x 10^ K. Therefore this absorption could arise in the warm/hot, 
shock-heated phase predicted by cosmological simulations. This is an upper limit because the profile 
could be broadened by factors other than thermal motions such as turbulence, multiple blended 
components, or expansion of the Universe. If the broadening is predominantly due to cosmological 
expansion, then the implied pathlength I through the absorber is given by Av « 15.0 km s^^ (//200 
kpc)/i75. Taking Av sa FWHM — 142 km s^^ for the broad H I component, we obtain I ^ 1.9 h'^^ 
Mpc in this case. 

4. The O VI absorption is not aligned with the strongest portion of the H I profile (i.e., the saturated 
component at u « km s"'^). Rather, the O VI is reasonably well-aligned with the broad H I 
component at w = 72 km s^^n The O VI line is also apparently broad; the profile fit shown in Figure 



Oegerle et al. (2000) suggested that the H I Ly/3 and O VI centroids were offset but could not be sure of the magnitude of 
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Fig. 2. — Continuum-normalized absorption profile of the O VI A1032 line at ^abs = 0.1212 (upper histogram), 
and the spectral region where the corresponding O VI A1038 line is expected (lower histogram) , both plotted 
versus rest-frame velocity where v = km s~^ at ^;abs = 0.1212. Both profiles are derived from the coadded 

data from different channels and observation dates (see text). The O VI A1038 line is strongly blended with 
a Ly(5 line at Zsths 

= 0.225. Nevertheless, there is evidence that the O VI A1038 line is present. The dashed 
line overplotted on the A1032 profile shows the result of a single-component fit to the 1032 A line only. 
The velocity, Doppler parameter, and column density from this fit were used to predict the strength of the 
corresponding (blended) A1038 transition, and the predicted 1038 A line is overplotted with a dashed line on 
the lower profile. An absorption feature with the predicted velocity, stength, and width is readily apparent. 
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Fig. 3. — Cluster of Lya absorption lines near the redshift of the O VI absorber at Zabs — 0.1212. The 
relative (unnormalized) flux is plotted versus rest-frame velocity where w = km s^^ at Zabs = 0.1212, and 
the continuum placement used for the measurements in Tables |l| and || is shown with a solid gray line. Seven 
galaxies with measured redshifts are close to the line of sight at z « 0.1212; their velocities in this frame 
[Aw = cAz/(l + z)] are indicated with tick marks at the top of the panel, and the numbers above each tick 
indicate the projected distance of the galaxy from the sight line. The uncertainties in the galaxy redshifts 
are estimated to be 50-100 km s^^. 
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has b = 59_2o s ^ . This would imply a high temperature if entirely due to thermal motions and 
may be evidence of another source of broadening (see below) . The upper limit on the temperature of 
the O VI gas for 6 = 59 km s'^ is T < 3.3 x 10^ K. 

5. As noted above, no metals apart from O VI are significantly detected despite coverage of very strong 
resonance transitions of abundant elements (e.g., C II A1334, C III A977, and Si III A1206). There is 
a marginal feature near the expected wavelength of C IV A1548 (see Figure but its significance is 
less than 3(t and no feature is evident at the expected wavelength of C IV A1550. Consequently, we 
use the 4cr upper limit on C IV for all purposes of this paper, and we obtain A^(0 VI)/A^(C IV) > 3.0. 

4.1. Collisional Ionization 

Can collisional ionization produce the observed properties of this O VI absorber? In coUisionally 
ionized gas in equilibrium (we briefly discuss non-equilibrium models at the end of §4.2), the O VI ionization 
fraction is maximized at T « 300,000 K (Sutherland & Dopita 1993). The H I ion fraction is small at 
such temperatures, but nevertheless, by virtue of the great abundance of hydrogen, H I Lya absorption 
associated with hot O VI should be detectable. Therefore H I absorption lines which are rather broad due 
to thermal motions are expected to arise in the warm/hot IGM discussed by Cen & Ostriker (1999a), and 
the broad H I and O VI lines shown in Figure |l] seem to provide promising evidence for the warm/hot gas. 
Does this hold up under scrutiny? The temperature implied by the apparent breadth of the H I component 
aligned with the O VI is about right for collisional ionization. The lower limit on the O VI/C IV ratio 
provides a lower limit on the temperature: assuming the gas is in equilibrium, we find from the calculations 
of Sutherland & Dopita (1993) that iV(0 VI)/7V(C IV) > 3.0 requires T > lO^'^ K. Combined with the 
temperature upper limit from the broad H I component, we have lO^ "^ <T< 10^-^ K. 

As noted above, the single-component O VI fit shown in Figure ^ implies T < 3.3 x 10^ K. While it is 
possible to detect trace O VI absorption in gas with T ^ 10^ K, the corresponding H I absorption would be 
quite broad and inconsistent with the observed Lya profile. It seems more likely that the O VI is at least 
partially broadened by non-thermal motions or multiple components. This non-thermal broadening would 
affect associated H I as well but would not preclude substantial thermal broadening in the H I profile. For 
example, if we assume that the broad H I and O VI absorption lines arise in the same gas and that the 
non-thermal motions can be adequately described by a Gaussian profile, then we can express the 6— value 
as = -I- 2kT/m and solve for T and 6nt, the contribution from non-thermal motions. With 6(H I) = 85 
km s^^ and 6(0 VI) = 59 km s^^, we obtain in this way T = 2.4 x 10^ K and b^t = 57 km s~^. 

The ionization mechanism can also be usefully tested by considering the metallicity required by the 
model — a model which requires an excessively high metallicity may be unrealistic. We estimate the gas 
metallicity as follows. Using the usual logarithmic notation, [0/H] = log(0/H) — log(O/H)0, the oxygen 
abundance can be expressed as 

"0 
H 



the offset because of problems with the wavelength scale in the early FUSE data. These problems have been rectified, and the 
offset is now reliably determined. 
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where / is the ion fraction and (O/H)© is the solar oxygen abundance.|^ Since A^(0 VI) and A^(H I) are 
measured quantities and (0/H)q is fixed, [O/H] can be estimated by using an ionization model to estimate 
/(H I)//(0 VI). 

Assuming the broad H I component shown in Figure |^ originates in the same gas as the O VI 



les lllc 
tion y 



and H Then, taking 
the oxygen abundance 



absorption, we have log N{R I) = 13.78 and log iV(0 VI) 14.02 from Tables 
/(H I)//(0 VI) vs. T from Sutherland & Dopita (1993), we derive from equati 
as a function of T required to match the measured A^(H I) and A^(0 VI). This is shown in Figure IJ. The 
lower limit on the temperature from A^(0 'V1)/N{C IV) and the upper limit on T from the H I line width 
are also indicated in Figure ^; the hatched gray regions are not consistent with these constraints. From this 
figure, we see that within the temperature limits, the observed column densities can be reproduced in gas 
with -1.8 < [O/H] < -0.6, i.e., 0.02 - 0.25 solar metallicity. 

This allowed range of metallicities is in reasonable agreement with theoretical expectations for 
intergalactic gas. For example, the study of IGM metal enrichment by Cen & Ostriker (1999b) indicates 
that this range of metallicity could arise in regions with a wide range of overdensity 5 = p/ < p > (see the 
z = curve in their Figure 2), and Z ^ 0.1 Z© would not be surprising given the H I column density of the 
absorber studied here. The models of Aguirre et al. (2001), which use various semi-analytic prescriptions to 
explore how metals might be distributed within a previously computed hydrodynamic simulation, predict 
somewhat lower mean metallicities at z = 0. The mean metallicity is somewhat less than 0.02 solar for 
S < 100 in their supernova-driven wind models, for example. However, their calculation only includes 
enrichment from galaxies with baryonic mass > 3 x 10^" M©, and inclusion of smaller galaxies could increase 
the metallicity by a factor of a few. Furthermore, since this is a mean metallicity, there are certainly higher 
metallicity regions in the models, and some of these are likely to be consistent with the requirements derived 
above for the Zabs = 0.1212 absorber. Of course, it is possible that this absorber does not originate in a 
particularly low overdensity region - the presence of several galaxies close to the absorber redshift suggests 
that this absorber may be due to gas within a galaxy group, or at least a somewhat overdense filament. 
However, the metallicity range shown in Figure |4| is reasonable in this scenario as well: X-ray observations 
(e.g., Davis, Mulchaey, & Mushotzky 1999; Hwang et al. 1999) of galaxy groups indicate that the intragroup 
medium metallicity is consistent with the upper portion of the range derived for this absorber. 

The H I and O VI absorption line kinematics are intriguing in the context of coUisional ionization. The 
juxtaposition of a relatively narrow H I component at w km s^^ with a broad H I + O VI component 
separated by ~75 km s~^ is suggestive of shock-heating in two colliding clouds, with the O VI and broad H I 
arising in the post-shock gas. If we estimate the shock velocity by scaling the apparent velocity separation 
of the components by \/3 (to account for projection effects), we estimate that two clouds colliding at a 
velocity = %/3 x 75 km s"^ = 130 km s'^ would be shock-heated to T w 2.3 x 10^ K, assuming a standard 
shockp] with specific heat ratio 7 = 5/3. This temperature, shown with a vertical dotted line in Figure ^, 
is consistent with the temperature constraints derived above.P| 



^"Throughout this paper we adopt the solar abundances reported by Grevesse & Anders (1989) and Grevesse & Noels (1993). 
These reference abundances have been widely used in recent literature. However, use of the recent revisions of the solar oxygen 
abundance (Holweger 2001; AUende Prieto et al. 2001; see also Sofia & Meyer 2001) increases [O/H] by ~ 0.1 - 0.2dex. 

^^For a standard shock with 7 = 5/3, the mean post-shock temperature is given by Tg = 1.38 X 10^(vs/100 km s~^)'^ K (see, 
e. g., McKee & HoUenbach 1980; Draine & McKee 1993). 

^^If the scale factor to correct for projection effects is less than Vs, then Ts could drop below the lower limit on the gas 
temperature set by N{0 Yl)/N{C IV). However, the component of the shock velocity which is projected onto the line-of-sight 
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log T (K) 

Fig. 4. — The logarithmic oxygen abundance [0/H] required to produce the observed O VI and H I column 
densities in the broad component of the absorber at Zabs = 0.1212, assuming the gas is coUisionally ionized 
and in equilibrium, as a function of temperature (solid line) . This follows from the O VI and H I ionization 
fractions vs. T (from Sutherland & Dopita 1993) and equation |l|. The gray regions of the plot indicate regions 
which are not allowed by other constraints; the lower limit on T is set by the lower limit on A^(0 VI) /N{C IV), 
and the upper limit on T is set by the 5— value of the broad H I component associated with the O VI 
absorption. The vertical dotted line indicates the mean post-shock temperature for a standard shock with 
Vs = 130 km s^^. 
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Summarizing this section, we find that the observed properties of this absorber can be produced 
by coUisionally ionized gas in equilibrium with 5.3 < log T < 5.6 and —1.8 < [0/H] < —0.6. These 
properties are consistent with expectations for the IGM as well as the intragroup medium in galaxy groups. 
Equilibrium coUisional ionization is evidently a viable ionization mechanism for this absorber. 



4.2. Photoionization 

Is photoionization also viable for this absorber? To address this question we have constructed standard 
models with the photoionization code CLOUDY (Ferland et al. 1998). The models adopt the usual 
assumptions: the absorber has constant density and the geometry is plane-parallel, the gas is photoionized 
by the UV background from QSOs and AGNs at z w 0.12 as calculated by Haardt & Madau (1996) ,0 
and the mean intensity at 1 Rydberg is set to J^{LL) = 1 x 10"^'^ ergs s~^ cm~^ Hz~^ sr~^, a value in 
reasonable agreement with observations (e.g., Kulkarni & Fall 1993; Maloney 1993; Vogel et al. 1995; 
Donahue, Aldering, & Stocke 1995; Tumlinson et al. 1999; Dave & Tripp 2001) and theoretical predictions 
(e.g., Haardt & Madau 1996; Fardal, Giroux, & ShuU 1998; Dave et al. 1999; Shull et al. 1999b). With 
these assumptions, we varied the gas metallicity and ionization parameter U (= nj/nn = ionizing photon 
density/total hydrogen number density) to find the parameter ranges consistent with the observed O VI 
and H I column densities in the broad component and the upper limits in Table ^ Note that since 
J,y(LL) is fixed at an assumed value, varying U is tantamount to varying nu- Since the photoionization 
model computes the H ion fraction, once U is constrained, the total H column density and thickness (or 
pathlength) of the absorber are also constrained. 

The results of the photoionization models are shown in Figure ^ where we have used the /(O VI) and 
/(H I) from the photoionization code to derive [0/H] vs. log U using equation |l|. From the measurements 
available in Table 0, we find that the lower limit on A''(0 Y1)/N{C IV) is the most useful for constraining 
the ionization parameter; from this ratio we derive log U > —1.20. A rough upper limit on U is provided 
by the required pathlength through the absorber. If the pathlength is too large, the lines would be broader 
than the observed lines due to cosmological expansion. Using the rough upper limit on / derived in §4, we 
find log U <^ —0.6 . As before, the gray hatching in Figure ^ indicates regions which are not consistent with 
these constraints. From this figure, we see that the observational constraints are satisfied if —1.1 < [0/H] 
< —0.3 and 85 < I 1900 kpc. To derive these metallicities, we have assumed log A^(H I) = 13.78 from the 
broad component fit shown in Figure 1. In the photoionized scenario, the H I Lya absorption line arising 
in the O VI gas could be substantially narrower, and this would lead to a lower H I column density. This, 
in turn, would require an even higher oxygen abundance to produce the observed A^(0 VI) [see eqn. 1]. 
However, log A^(H I) — 13.78 is likely a reasonable upper limit in this case, and therefore these models 
provide lower limits on the metallicity if the gas is photoionized. 

The physical properties implied by the photoionization models are reasonable. For example, Schaye 
(2001) has analytically derived expressions for the density and size of Lya absorbers, assuming only that 
the gas is photoionized and in local hydrostatic equilibrium. Using his equations (8) and (12), and assuming 
the gas temperature T ~ lO'^K, the H I photoionization rate F 10"^^ '^ s~^, and the fraction of the mass in 



could also be larger than 75 km s ^, which would offset this effect. The point here is that the observed absorption lines have 
approximately the right velocities to produce O VI in a shock. 

^^We have also computed models using the UV background calculated by Shull et al. (1999b), and we obtain very similar 
results. 
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Fig. 5. — The logarithmic oxygen abundance [0/H] (sohd hne) required to produce the observed O VI 
and H I column densities in the broad component of the absorber at Zabs = 0.1212, assuming the gas is 
photoionized by the UV background from QSOs, as a function of the ionization parameter (bottom axis) 
and the total H density (top axis) for Jj.(H I) = 1 x lO'^^ ergs s^^ cm'^ Hz"! sr'^ The the O VI and H I 
ionization fractions vs. log U were calculated with CLOUDY (see text). The dashed line shows the required 
pathlength through the absorber vs. log U using the scale on the right axis. 



- 18 - 



gas /g ^ 0.2, we obtain log ?ih ^ -5.6 and I ^ 300 kpc for fog A^(H I) < 13.78. These properties are entirely 
consistent with the properties inferred from the CLOUDY models (see Figure [s]). Perhaps the most serious 
concern about these photoionization models is that the IGM must be substantially and uniformly enriched 
with metals over a long pathlength. This problem is less severe at higher metallicities because the required 
pathlengths are smaller. However, the minimum absorber thickness (set by the lower limit on U from 
O VI/C IV) is 85 kpc for [0/H] — —0.3. Some QSO absorbers have been observed which apparently have 
very low metallicities (e.g., ShuU et al. 1999a), so the IGM has not been generally enriched to this level, but 
the system of this paper could arise in a more metal-rich region. As noted above, the mean metallicities 
predicted by simulations (e.g., Gen & Ostriker 1999b; Aguirre et al. 2001) are only marginally consistent 
with such a high metallicity, but again this could be a particularly high Z region. Gen & Bryan (2001) have 
argued that most of the metals in the diffuse, low-density intergalactic gas are provided by small galaxies 
{M < IO^Mq) and are injected into the IGM at z > 4. Because larger galaxies which form later do not 
widely disperse their metals, Cen & Bryan expect that low-z, low column density Lya clouds will have low 
metallicities. A similar picture has been favored by Heckman et al. (2000) based on the observed properties 
of supernova-driven "superwinds" from low— 2; star burst galaxies and high— 2: Lyman-break galaxies (Pettini 
et al. 2001). Heckman et al. argue that metals injected into the IGM via superwinds, primarily from the 
lower-mass galaxies, could increase the IGM metallicity to ~l/6 solar metallicity in the intracluster medium 
as well as the general IGM. However, in most supernova-wind enrichment scenarios, the O VI gas should be 
predominantly coUisionally ionized. Radiation-pressure driven outflows of dust from galaxies (e.g., Aguirre 
et al. 2001) could transport heavy elements more gently. 

It is interesting to note that Savage et al. (2001) have identified an intervening O VI absorber at Zabs 
= 0.06807 in the spectrum of PG0953-I-415 which is apparently well-explained by photoionization. In this 
case, H I, C III, G IV, and N V are also detected and are well-aligned with the O VI, and all of the column 
densities are reasonably matched by a photoionization model, albeit with a somewhat uncomfortably high 
metallicity /pathlength combination. Perhaps more importantly, the H I Lya and Ly/3 lines associated with 
the O VI are narrow and appear to be composed of only one component, and the upper limit on T from 
the width of the H I lines is T < 4.1 x lO'' K. This precludes collisional ionization, at least in equilibrium, 
unless the H I arises in a different phase from the O VI. A similar situation is found for the O VI system 
studied by Tripp & Savage (2000), although in this case it is possible to hide a broad component in the 
complicated Lya profile, assuming the absorber is a multiphase medium (see their Figure 6). While these 
examples would seem to strongly favor photoionization, it is also possible that they represent cases where 
the gas was shock-heated to ^ 10^ K and then cooled more rapidly than it could recombine, as discussed 
by various authors (e.g., Edgar & Chevalier 1986). Table 4 in Tripp & Savage (2000) summarizes high ion 
ratios predicted by four non-equilibrium collisional ionization models. Comparison of the high ion column 
densities measured in this paper to that table shows that most non-equilibrium collisional ionization models 
satisfy the observational constraints. Of course, in this hypothesis, the substantial breadth of the H I profile 
would be unexpected, and some other broadening mechanism would need to be invoked. This could simply 
be the kinematics of the gas, e.g., in a galactic wind. 



4.3. X-ray Absorption Lines 

In principle, one promising means to break the collisional ionization/photoionization degeneracy is 
to search for X-ray absorption lines, primarily the resonance transitions of O VII and O VIII at 574 eV 
and 654 eV, respectively (21.6 and 19.0 A). Several groups have discussed various aspects of the use of 
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these lines and their detection with current X-ray facihties (Shapiro & Bahcall 1980; Aldcroft et al. (1994); 
Hellsten, Gnedin, & Miralda-Escude 1998; Perna & Loeb 1998; Fang & Canizares 2000; Fang et al. 2001). 
A complete summary of this literature is beyond the scope of this paper, but we note that the first attempt 
to detect these absorption lines in the diffuse IGM (with the Chandra High Energy Transmission Grating 
Spectrometer) has yielded upper limits of A^(0 VIII) < 10^^ cm~^ (Fang et al. 2001). 

Given this rough guideline for detection capabilities with current X-ray telescopes, can we hope to 
determine the ionization mechanism of these O VI absorbers by searching for the corresponding X-ray 
absorption lines at the same redshifts? To explore this question, we show the O VI/0 VII, O VI/0 VIII, and 
O VII/0 VIII column density ratios predicted for collisionally ionized gas (in equlibrium) and photoionized 
gas in Figures ^ and ||b, respectively, for the expected ranges of T and U. The ratios for collisionally 
ionized gas are based on the calculations of Sutherland & Dopita (1993). The ratios for photoionized gas 
are from the CLOUDY models described in §4.2. Equilibrium collisional ionization and photoionization 
clearly have different signatures in these ratios. For example, the O VI/0 VII ratio has a much broader 
range in collisionally ionized gas than in photoionized gas: given the temperature and ionization parameter 
constraints derived in §4 for the particular absorber of this paper, the models predict — 1 ^ log (O VI/0 VII) 
< 1 for collisionally ionized gas and log (O VI/0 VII) w for photoionized gas. Furthermore, if log 
(O VI/0 VII) « 0, then the O VIII line can distinguish between the two models because in this situation 
the O VI/0 VIII ratio is much larger in the collisionally ionized model than the photoionized model.[^ 
While Figure ^ demonstrates the potential of X-ray observations, it unfortunately also shows that for the 
O VI system studied here, the O VII and O VIII lines are predicted to have column densities well below 10^^ 
cm^^ and therefore may be difficult or impossible to detect with current instruments. X-ray observations of 
other absorbers with larger O VI column densities may be more illuminating. However, if the O VI column 
density is too much larger, then the O VI lines will saturate making iV(0 VI) difficult to reliably measure. 
It seems that it will be difficult to study O VI, O VII, and O VIII absorption arising from the same gas. If 
the higher oxygen ions were detected at the redshift of one of the H1821-I-643 O VI systems, this would 
likely indicate the presence of multiphase gas with the O VII and O VIII revealing a more highly ionized 
phase containing little O VI. 

This discussion has several caveats. First, we have assumed that either collisional ionization or 
photoionization dominates. In fact, both ionization mechanisms may contribute. Hellsten et al. (1998) 
argue that photoionization by the X-ray background dominates in the production of O VII and O VIII, 
even if the gas has a temperature of 10^ — 10*" K, but this depends on the density of the gas. The O VI 
might also be produced by the combined effects of collisional and photoionization. Second, the gas is 
assumed to be in ionization equilibrium. We have already noted that this assumption might not hold for 
collisionally ionized gas. Photoionized O VI is likely to be in equilibrium, but O VII and O VIII produced 
by photoionization might be out of equilibrium because the ionization time scales are comparable to the 
Hubble time (see, e.g., Hellsten et al. 1998). Again, it would be valuable to evaulate the ratios predicted for 
these ions in non-equilibrium models as well. Finally, the intensity of the X-ray background is uncertain. A 
full assessment of these issues is beyond the scope of this paper. 



Along these lines we note that in general, the O VII and O VIII lines without the O VI constraint may not be able to 
distinguish between collisional ionization and photoionization because the dotted lines in Figure ^ are degenerate, and the 
resolution of current X-ray spectrometers is not sufficient to constrain the temperature of the gas based on line widths. 
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Fig. 6. — Logarithmic column density ratios predicted for O VI/0 VII (solid line), O VI/0 VIII (dashed 
line), and O VII/0 VIII (dotted line) for (a) coUisionally ionized gas in equilibrium, and (b) gas photoionized 
by the UV background from QSOs. The ratios in the coUisionally ionized case are plotted vs. log T, while 
the photoionized gas ratios are shown vs. log of the ionization parameter. 
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5. Summary and Discussion 

High-resolution STIS and FUSE UV spectra of H1821+643 reveal that the absorption system at 
^abs — 0.1212 has a number of interesting properties including (1) detection of O VI and H I absorption 
without significant detection of any other species, (2) complex H I absorption including an apparently 
broad component with 6 « 85 km s~^ that is aligned with the O VI absorption, and (3) association with 
one or more galaxies close to the sight line. Using constraints from these spectra, we have examined the 
ionization mechanism in this absorber. We find that equilibrium collisional ionization is viable with 5.3 < 
log T < 5.6 and —1.8 < [0/H] < —0.6. However, the absorption line properties can also be explained by 
photoionization if —1.1 < [0/H] < —0.3 and 85 < ^ 1900 kpc. In addition, it is possible that the gas is 
collisionally ionized but is not in equilibrium (e.g., Edgar & Chevalier 1986). 

There are a number of possible sites in which O VI absorption lines could arise in the spectra of 
low— z QSOs. In some cases, O VI absorption likely occurs in the immediate environment of the QSO itself 
(e.g.. Savage et al. 1998; Papovich et al. 2000; Ganguly et al. 2001, and references therein). The large 
displacement from the QSO redshift and the close association with galaxies near the line of sight indicates 
that the absorber studied in this paper is not in this class. Rather, it is an intervening absorber. However, 
the detailed nature of the intervening absorption is still unclear - it could occur in the bound ISM of an 
individual galaxy, the escaping gas in a galactic wind, the intragroup medium of a galaxy group, or a more 
remote region of a large scale structure or the diffuse IGM. 

The galactic wind hypothesis is appealing because O VI absorption has been detected in such an outflow 
from the nearby starburst galaxy NGC1705 (Heckman et al. 2001). Gas in the temperature range needed to 
produce O VI by collisional ionization is predicted in hydrodynamic simulations of starburst galaxy winds 
(e.g., Mac Low & Ferrara 1999; Strickland & Stevens 2000). However, unlike the O VI absorber studied in 
this paper, the outflow from NGC1705 shows strong absorption lines due to low and intermediate ionization 
stages as well as O VI. This absorption by low ions is consistent with the hydrodynamic models in which 
the O Vl-bearing gas arises in interfaces between cool, denser gas and much hotter gas (see Figure 4c in 
Heckman et al. 2001). In this model, low ion absorption occurs in the cool gas, O VI exists in the interface, 
and the hotter gas produces X-ray emission. The NGC1705 observations probe the ISM and outflow close 
to the central starburst and its host galaxy, but the Zabs — 0.1212 O VI system is at a substantially larger 
projected distance from any galaxy known to be at the absorber redshift. One could speculate that the 
character of absorption from a wind is different at larger distances away from the star cluster, but these 
winds have not been observationally or theoretically studied at such large impact parameters. 

The hypothesis that the absorption occurs in the ISM of a galaxy like the Milky Way suffers similar 
criticisms. While O VI is commonly detected in the disk and halo of the Milky Way (Jenkins 1978; Savage 
et al. 2000), it would typically be accompanied by absorption from other heavy elements. Galactic high 
velocity clouds, including clouds stripped out of satellite galaxies, also frequently show associated O VI 
(Sembach et al. 2000), and high velocity clouds have approximately the same metallicity as the absorber 
studied in this paper (Wakker et al. 1999; Gibson et al. 2000; Richter et al. 2001). However, a sight line 
through such an entity would have a much higher H I column density and again would be expected to be 
detected in other ions.^ 

Perhaps the most appealing hypothesis is that the Zabs — 0.1212 absorber is due to gas between galaxies 



'^^The particular absorber studied here is only detected in O VI and H I, but other intervening O VI systems have been 
detected in other species (Tripp et al. 2000; Chen & Prochaska 2000; Savage et al. 2001). 
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in a group or unvirialized filamentary structure. Diffuse X-ray emission indicating the presence of a liot 
intragroup medium has been detected from many galaxy groups (Mulchaey 2000 and references therein), 
but mainly groups that are elliptical-rich (Mulchaey et al. 1996a; Zabludoff & Mulchaey 1998). Mulchaey 
et al. (1996b) have suggested that spiral-rich galaxy groups also contain hot intragroup gas, but the gas 
is not quite hot enough to be detected in X-ray emission (T < 4 x 10^ K). They predict that the O VI 
absorption lines should be detectable in such an intragroup medium. It is interesting that the Zabs absorber 
fits this hypothesis in many respects: (1) there are several galaxies near the sight line at this redshift (see 
Figure H), (2) at least one of the galaxies is a spiral/actively star-forming galaxy, and (3) the breadth of the 
O VI absorption is considerably larger than expected from thermal motions (§4), which is consistent with 
additional broadening due to the velocity dispersion of the group.^ However, this may not be a bound 
group. Given our currently limited information, it is difficult to derive additional information about the 
collection of galaxies near this O VI absorber, but we do note that the spatial extent of the ensemble appears 
to be larger than that of a typical poor group (c.f., Zabludoff & Mulchaey 1998). It would be valuable to 
obtain additional galaxy redshift measurements and high-resolution images for morphological classification. 

In §1 we outlined several areas in which additional effort is needed to elucidate the nature of the O VI 
systems and their contribution to the baryon budget. The main goal of this paper has been to scrutinize the 
ionization mechanism to test the idea that a substantial quantity of shock-heated hot gas is present in the 
IGM at the current epoch. We find that the O VI absorber at Zabs = 0.1212 in the spectrum of H1821+643 
has properties consistent with coUisionally ionized hot gas such as a broad H I component aligned with 
the O VI and a high O VI/C IV ratio. This is encouraging evidence in support of the shock-heated IGM 
hypothesis. However, we cannot rule out photoionization or a combination of coUisional and photoionization 
processes. The degree of IGM metal enrichment required by photoionization exceeds the predictions of some 
models, but given our limited understanding of how metals are transported out of galaxies and the variety 
of possibilities for the absorption site, this is not sufhcient grounds to dismiss photoionization. Analyses of 
other low— z O VI systems (e.g., Tripp & Savage 2000) have reached similar conclusions regarding coUisional 
vs. photoionization. 

Given the general difhculty encountered in attempts to definitively identify the ionization mechanism 
in individual systems, it will be important carry out statistical analyses of O VI absorber samples including 
comparisons to various models. For example. Gen et al. (2001) and Fang & Bryan (2001) have recently 
predicted the properties of O VI absorbers based on cosmological simulations. The absorber statistics 
predicted by these models [e.g., dN/dz and r2b(0 VI)] appear to be in reasonable agreement with current 
observations. In these simulations, O VI absorption arises in both coUisionally ionized and photoionized gas. 
However, the photoionized systems are narrower and have lower equivalent widths than the coUisionally 
ionized absorbers. Tripp (2002) has recently presented the 6— values and column densities of twenty low— z 
O VI lines observed with STIS and FUSE (see his Figure 5). It is interesting to note that while there are a 
few apparently narrow O VI lines in the observations compiled by Tripp (2002), the majority of the O VI 
lines in that sample have 6— values consistent with an origin in hot gas. The &— value of an O VI line in gas 
at T ~ 300, 000 K is ~18 km s~^, and the median &— value of the observed sample is 22 km s""'^. In addition, 
the apparently narrow O VI lines tend to be weaker systems as well, as predicted by the simulations. These 
preliminary results appear to be fully consistent with the cosmological simulations. We look forward to 
similar analyses using larger observational samples and a variety of models. 



The typical velocity dispersion of a spiral-rich group is ~ 100 km (Mulchaey et al. 1996b). 
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